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Abstract 

I abstract from a recent publication [T] the motivations for, analysis in and 
conclusions of a study of the ultraviolet and infrared momentum regulators in- 
duced by the necessary truncation of the model spaces formed by a variational 
trial wave function. This trial function is built systematically from a complete 
set of many-body basis states based upon three-dimensional harmonic oscillator 
(HO) functions. Each model space is defined by a truncation of the expansion 
characterized by a counting number [M) and by the intrinsic scale (fio;) of the 
HO basis. Extending both the uv cutoff to infinity and the ir cutoff to zero is 
prescribed for a converged calculation. In [I] we established practical procedures 
which utilize these regulators to obtain the extrapolated result from sequences 
of calculations with model spaces. Finally, I update this subject by mentioning 
recent work on our extrapolation prescriptions which have appeared since the 
submission of [l]. The numerical example chosen for this contribution consists 
of calculations of the ground state energy of the triton with the "bare" and 
"soft" Idaho N'^LO nucleon-nucleon (NN) interaction. 

Keywords: no-core shell model; convergence of expansion in harmonic oscilla- 
tor functions; ultraviolet regulator, infrared regulator 

1 Introduction 

The advent of giant nuclear shell-model codes based upon the three-dimensional har- 
monic oscillator (HO) in the 1970s coincided with the advent of a program to use 
the HO eigenf unctions as a basis of a finite linear expansion to make a straightfor- 
ward variational calculation of the properties of light nuclei [5]. At the same time 
theorems based upon functional analysis established the asymptotic convergence rate 
of these latter calculations as a function of the counting number (call it J\f) which 
characterizes the size of the expansion basis (or model space) O H]. The conver- 
gence rates of these theorems (inverse power laws in Af for "non smooth" potentials 
with strong short range correlations and exponential in A/" for "smooth" potentials 
such as gaussians) were demonstrated numerically in |3] for the HO expansion and in 
[5] for the parallel expansion in hyperspherical harmonics. These convergence theo- 
rems seem to be known in the hyperspherical harmonic community and are effectively 
demonstrated in the calculation of the properties of few-nucleon systems "from first 
principles" ; that is, solving the many-body Schrodinger equation with a Hamiltonian 
containing nucleon-nucleon interactions fitted to scattering data and to properties of 
the deuteron. The convergence rates of variational calculations using the HO basis 
have been periodically rediscovered empirically by those who, in the present day, have 
adapted "giant shell- model codes" or written new codes to perform "a& initio'' "no- 
core shell model" (NCSM) calculations of s- and p-shell nuclei. I have never seen a 
reference to the functional analysis theorems regarding these convergence rates in the 
NCSM papers. However, the HO expansion basis has an intrinsic scale parameter huj 
which does not naturally fit into an extrapolation scheme based upon Af as discussed 
by [3 m E] . Indeed the model spaces of these NCSM approaches are characterized 
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by the ordered pair {Af,huj). Here the basis truncation parameter Af and the HO 
energy parameter hw are variational parameters [2l 13 [9] . It is the purpose of this 
contribution to summarize the properties of another ordered pair which perhaps more 
physically describes the nature of the model spaces and provides extrapolation tools 
which use Af and hui on an equal footing [T]. This is the pair of ultraviolet (uv) and 
infrared (ir) cutoffs (each a function of both Af and Hlo) induced by the truncation. 
They were first introduced to the NCSM in [10] in the context of an effective field 
theory (EFT) approach (for a recent review of this program see [TT]) . These cutoffs 
or regulators can usefully be employed in novel extrapolation schemes [I] which are 
a natural outgrowth of those introduced in the 1970s, rediscovered by the NCSM 
community, and in current use. 

The variational approach alluded to above generates a trial wave function in a 
completely systematic manner without regard for the details of the Hamiltonian under 
consideration other than the implementation of exact symmetries. The goal, then, 
is to define a complete set of states for a few-body system and to construct and 
diagonalize the Hamiltonian matrix in a truncated basis of these states. The result 
of the diagonalization is an upper bound to the exact eigenvalue of the complete 
set. With this method a reliable estimate of the accuracy attained can be made 
with the variational upper bound [3] provided that the trial function is constructed 
using the terms of a systematic expansion set and convergence of the diagonalization 
result (such as a ground-state energy) is observed as the basis is increased. The 
algebra appropriate to generating and using trial wave functions, based on three 
dimensional HO eigenfunctions, has been given by Moshinsky [2] and others [12] ■ The 
trial functions take the form of a finite linear expansion in a set of known functions 

V 

where a^^^ are the parameters to be varied and hy are many-body states based on 
a summation over products of HO functions. The advantage of a HO basis is that 
it is relatively straightforward to construct a complete set of few-body functions of 
appropriate angular momentum and symmetry; examples are given in [121 113| . The 
trial function must have a definite symmetry reflecting the composition of the bound 
state: fermions or bosons. This trial function '^t must be quadratically integrable and 
the expectation value of the Hamiltonian must be finite. The expansion coefficients 
(known as generalized Fourier coefficients in the mathematical literature) depend 
on the upper limit (such as an defined in terms of total oscillator quanta) and 
are obtained by minimizing the expectation value of the Hamiltonian in this basis. 
Treating the coefficients al/^-* as variational parameters in the Rayleigh quotient , 
one performs the variation by diagonalizing the many-body Hamiltonian in this basis. 
This is an eigenvalue problem so the minimum with respect to the vector of expansion 
coefficients always exists and one obtains a bound on the lowest eigenvalue (and indeed 
on the higher eigenvalues representing the excited states [H]). The basis functions 
can also depend upon a parameter (such as the harmonic oscillator energy hio which 
sets a scale) that then becomes a non-linear variational parameter additional to the 
linear expansion coefficients. 

One can view a shell-model calculation as a variational calculation, and thus ex- 
panding the configuration space merely serves to improve the trial wave function |16j . 
The traditional shell-model calculation involves trial variational wave functions which 
are linear combinations of Slater determinants. Each Slater determinant corresponds 
to a configuration of A fermions distributed over A single-particle states. If we take 
any complete set of orthonormal single-particle wave functions and consider all pos- 
sible ^-particle Slater determinants that can be formed from them, then these wave 
functions form a complete orthonormal set of wave functions spanning the A-particle 
Hilbert space. The Slater determinant basis of HO single-particle wave function is 
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often defined in the "m-scheme". That is, the single-particle states are labelled by 
the quantum numbers n, /, j, and rrij, where n and I are the radial and orbital HO 
quantum numbers, j is the total single-particle spin, and mj its projection along the 
z-axis. The many-body basis states have well-defined total spin projection, which 
is simply the sum of mj of the single-particle states Mj — ^ rrij , hence the name 
"m-scheme" . The many-body basis states are limited only by the imposed symme- 
tries — parity, charge and total angular momentum projection (M) , as well as by Af. 
However, in general the many-body basis states do not have a well defined total J. 
This scheme is simple to implement and in two calculations (for positive and negative 
parity) one gets the complete low-lying spectrum, including the ground state, even 
though the spins of the low-lying states are not specified in the trial wave function. 
The truncation by Af results in finite matrices to be diagonalized, but they are much 
larger than the matrices of the Moshinsky program which expects the properties of the 
trial wave function ( JT basis in shell model language) to be known. However, because 
these shell model wave functions do span the space, an expansion in such "m-scheme" 
Slater determinants is, in principle, also capable of giving an exact representation of 
the eigenfunctions of the Hamiltonian. 

These early ab initio calculations, both of the "no-core" shell model in which all 
nucleons are active [1^] and of the Moshinsky program [T71 [TB] attempted to overcome 
the challenges posed by "non-smooth" two-body potentials by including Jastrow type 
two-body correlations in the trial wave function. Nowadays, the NN potentials are 
tamed by unitary transformations within the model space |19j or in free space by 
the similarity renormalization group evolution j20j . In both cases, this procedure 
generates effective many-body interations in the new Hamiltonian. Neglecting these 
destroys the variational aspect of the calculation (and the physics contained in the 
calculation, of course). We retain the variational nature of our NCSM investigation 
by choosing a realistic smooth nucleon-nucleon interaction Idaho N'^LO [21] which has 
been used previously without renormalization for light nuclei(A < 6) [71. This poten- 
tial is inspired by chiral perturbation theory and fits the two body data quite well. It 
is composed of contact terms and irreducible pion-exchange expressions multiplied by 
a regulator function designed to smoothly cut off high-momentum components in ac- 
cordance with the low-momentum expansion idea of chiral perturbation theory. The 
version we use has the high- momentum cutoff of the regulator set at 500 MeV/c. The 
Idaho N'^LO potential is a rather soft one, with heavily reduced high-momentum com- 
ponents as compared to earlier realistic NN potentials with a strongly repulsive core. 
Alternatively, in coordinate space, the Yukawa singularity at the origin is regulated 
away so that this potential would be considered "smooth" by Delves and Schneider 
and the convergence in Af would be expected to be exponential [SI E] . Even without 
the construction of an effective interaction, convergence with the Idaho N^LO NN 
potential is exponential in Af, as numerous studies have shown [71 [5^ . 

With the HO basis in the nuclear structure problem, convergence has been dis- 
cussed, in practice, with an emphasis on obtaining those parameters which appear 
linearly in the trial function (i.e. convergence with Af). Sometimes for each Af the 
non-linear parameter huj is varied to obtain the minimal energy [7, [22j for a fixed Af 
and then the convergence with Af is examined. Sometimes Huj is simply fixed at a 
value which gives the fastest convergence in A/" [13] . Other extrapolation schemes have 
been proposed and used [8]. In all of these schemes, in my opinion, the extrapolation 
to an infinite basis is effected with the main role played by Af and a secondary role 
played by huj. The scheme proposed in 1 gives Af and huj equal roles by employing 
uv and ir cutoffs which which must be taken to infinity and to zero, respectively to 
achieve a converged result (see Figure 1). 

In section 2 we briefly describe expansion schemes in HO functions. None of the 
discussion in section 2 is new, but it paves the way for section 3 in which we suggest 
a convergence analysis based upon the uv and ir cutoffs induced by the truncation 



4 



Sidney A Coon 



A 



Ultraviolet cutoff 




A 



Excluded due to finite N '^"^-.^ 
(due to lack of computer 
power) 



/ 



s 

c 



/ 




Low-energy 
physics content of 
the nucleus 



o 



E 

3 
C 

a» 
E 




Content of the nucleus obtained by 
"ab-initio" no-core shell model methods 



Excluded due to finite N 
(due to lack of computer 
power) .y 



i 



V 



\ 



Infrared cutoff 




V 



Figure 1: (Color online) Schematic view of a finite model space (limited by the basis 
truncation parameter N as described in the text), in which the uv and ir momentum 
cutoffs are arbitrary. To reach the full many-body Hilbert space, symbolized by the 
complete oval, one needs to let the uv cutoff -> oo and the ir cutoff — 0. 

of the model space. Section 4 is devoted to a sampling of tests and examples of this 
new convergence scheme; for a more extensive discussion with more examples please 
see Ref. [T]. 

2 Expansion in a finite basis of harmonic oscillator 
functions 

We briefly indicate the workings of the finite HO basis calculations performed and re- 
fer the reader to a comprehensive review article f5] on the no-core shell model (NCSM) 
for further details and references to the literature. A HO basis allows preservation of 
translational invariance of the nuclear self-bound system. Translational invariance is 
automatic if the radial HO wave function depends on relative, or Jacobi, coordinates 
as was done in Refs. [13l|T7l|T8l|22] . Antisymmetrization (or symmetrization for the 
a particle models of [T71[TH]) of the basis is necessary and described in Refs. [3] and 
[23] . Antisymmetrization in a Jacobi basis becomes analytically and computationally 
forbidding as the number of nucleons increases beyond four or five. For this rea- 
son these calculations are alternatively made with antisymmetrized wave functions 
constructed as Slater determinants of single-nucleon wave functions depending on 
single-nucleon coordinates. This choice loses translational invariance since, in effect, 
one has defined a point in space from which all single-particle coordinates are defined. 
Translational invariance is restored by choosing a particular truncation of the basis: a 
maximum of the sum of all HO excitations, i.e. ^f^ii'^ni + k) < Ntotmax, where n^, k 
are the HO quantum numbers corresponding to the harmonic oscillators associated 
with the single-nucleon coordinates and Ntotmax is an example of the generic Af of 
the Introduction. The gain of this choice is that one can use technology developed 
and/or adapted for NCSM, such as the shell model code ANTOINE [IS], the parallel- 
processor codes "Many-Fermion Dynamics — nuclear" (MFDn) [25] and the No-Core 
Shell Model Slater Determinant Code [57] • These codes set up the many-body basis 
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space, evaluate the many-body Hamiltonian matrix, obtain the low-lying eigenvalues 
and eigenvectors using the Lanczos algorithm, and evaluate a suite of expectation 
values using the eigenvectors. 

The eigenstates factorize as products of a wave function depending on relative 
coordinates and a wave function depending on the cm. coordinates. The precise 
method of achieving the factorization of the cm. and intrinsic components of the 
many-body wave function follows a standard approach, sometimes referred to as the 
"Lawson method" [21]. In this method, one selects the many-body basis space in 
the manner described above with Af = Ntotmax and adds a Lagrange multiplier term 
to the many-body Hamiltonian l5{Hc.m. — where Hem. is the HO Hamiltonian 

for the cm. motion. With /3 chosen positive (10 is a typical value), one separates 
the states of lowest cm. motion (05*1 ) from the states with excited cm. motion by 
a scale of order j3hio. The resulting low- lying states have wave functions that then 
have the desired factorized form. We checked, for the two cases ^ = 3 and A — A, 
that the codes manyeff [23] which use Jacobi coordinates and No-Core Shell Model 
Slater Determinant Code ^\ based upon single-nucleon coordinates gave the same 
eigenvalues for the same values of N = Ntotmax and hio, indicating that the Lawson 
method is satisfactory for the calculations in single-particle coordinates. 

Now we return to the truncation parameter N of the HO basis expansion of the 
many-body system. Usually, instead of truncating the sum of all HO excitations 
N — Ntotmax, one uses the the more familiar truncation parameter Nmax- Nmax is 
the maximum number of oscillator quanta shared by all nucleons above the lowest 
HO configuration for the chosen nucleus. One unit of oscillator quanta is one unit of 
the quantity (2n -|- I) where n is the principle quantum number and I is the angular 
quantum number. For A = 3,4 systems Nmax = Ntotmax- For the p-shell nuclei they 
differ, e.g. for "^Li, Nmax = Notmax-'^, and for ^^C, Nmax = NotmaxS. Later on we 
will want a truncation parameter which refers, not to the many-body system, but to 
the properties of the HO single-particle states. If the highest HO single-particle state 
of this lowest HO configuration has Nq HO quanta, then Nmax + Nq — N identifies 
the highest HO single-particle states that can be occupied within this many-body 
basis. Since Nmax is the maximum of the total HO quanta above the minimal HO 
configuration, we can have at most one nucleon in such a highest HO single-particle 
state with N quanta. Note that Nmax characterizes the many-body basis space, 
whereas is a label of the corresponding single particle space. Let us illustrate this 
distinction with two examples. ^He is an open shell nucleus with Nq — 1 since the 
valence neutron occupies the Op shell in the lowest many-body configuration. Thus 
if Nmax = 4 the single particle truncation TV is 5. On the other hand, the highest 
occupied orbital of the closed shell nucleus ^He has iVo = so that N — Nmax- 



3 Ultraviolet and infrared cutoffs induced by basis 
truncation 

We begin by thinking of the finite single-particle basis space defined by N and hu 
as a model space characterized by two momenta associated with the basis functions 
themselves. In the HO basis, we follow [TD] and define A = y/mpj{N + 3/2)huj as 
the momentum (in units of MeV/c) associated with the energy of the highest HO 
level. The nucleon mass is tojv = 938.92 MeV. To arrive at this definition one applies 
the virial theorem to this highest HO level to establish kinetic energy as one half the 
total energy {i.e., {N -I- 3/2)Sa; ) and solves the non-relativistic dispersion relation 
for A. This sets one of the two cutoffs for the model space of a calculation. En- 
ergy, momentum and length scales are related, according to Heisenberg's uncertainty 
principle. The higher the energy or momentum scale we may reach, the lower the 
length scale we may probe. Thus, the usual definition of an ultraviolet cutoff A in 
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the continuum has been extended to discrete HO states. It is then quite natural to 
interpret the behavior of the variational energy of the system with addition of more 
basis states as the behavior of this observable with the variation of the ultraviolet 
cutoff A. Above a certain value of A one expects this running of the observable with 
A to "start to behave" so that this behavior can be used to extrapolate to the exact 
answer. However, the truncation of the model space by Af implies a second cutoff, 
absent in free space; an infrared cutoff. Because the energy levels of a particle in a HO 
potential are quantized in units of Hlj, the minimum allowed momentum difference 
between single-particle orbitals is A = ^/mi^hhj and that has been taken to be an 
infrared cutoff [TU]. That is, there is a low-momentum cutoff A = fi/b corresponding 



characteristic length of the HO potential and basis functions). Note however that 
there is no external confining HO potential in place. Instead the only huj dependence 
is due to the scale parameter of the underlying HO basis. In [lOj the infiuence of the 
infrared cutoff is removed by extrapolating to the continuum limit, where hiu — )■ 
with — > oo so that A is fixed. Clearly, one cannot achieve both the ultraviolet limit 
and the infrared limit by taking tuo to zero in a fixed- A^model space as this procedure 
takes the ultraviolet cutoff to zero. 

The calculated energies of a many-body system in the truncated model space will 
differ from those calculated as the basis size increases without limit {N — oo). This 
is because the system is in effect confined within a finite (coordinate space) volume 
characterized by the finite value of h intrinsic to the HO basis. The "walls" of the 
volume confining the interacting system spread apart and the volume increases to the 
infinite limit as A — >■ and & — >■ oo with A held fixed. Thus it is as necessary to 
extrapolate the low momentum results obtained with a truncated basis with a given h 
or huj as it is to ensure that the ultraviolet cutoff is high enough for a converged result. 
These energy level shifts in a large enclosure have long been studied [55] ; most recently 
with the explicit EFT calculation of a triton in a cubic box allowing the edge lengths 
to become large (and the associated ir cutoff due to momentum quantization in the 
box going towards zero) [29]. There it was shown that as long as the infrared cutoff 
was small compared to the ultraviolet momentum cutoff appearing in the "pionless" 
EFT, the ultraviolet behavior of the triton amplitudes was unaffected by the finite 
volume. More importantly, from our point of view of desiring extrapolation guidance, 
this result means that calculations in a finite volume can confidently be applied to the 
infinite volume (or complete model space) limit. Similar conclusions can be drawn 
from the ongoing studies of systems of two and three nucleons trapped in a HO 
potential with interactions from pionless EFT combined with this definition of the 
infrared cutoff (A = y/m^fhuj); see the review 

Other studies define the ir cutoff as the infrared momentum which corresponds 
to the maximal radial extent needed to encompass the many-body system we are 
attempting to describe by the finite basis space (or model space). These studies find 
it natural to define the ir cutoff by Asc = a/ {'m^huj)/{N + 3/2) [23 133 • Note that A^c 
is the inverse of the root-mean-square (rms) radius of the highest single-particle state 
in the basis; (r^)^/^ = by^N + 3/2. We distinguish the two definitions by denoting 
the first (historically) definition by A and the second definition by Age because of its 
scaling properties demonstrated in the next Section. 

The extension in |10| of the continuum ultraviolet cutoff to the discrete (and trun- 
cated) HO basis with the definition A = yJmN{N + 3/2)hoj seems unexceptional. 
But, as always when one confidently makes such a statement, there are exceptions. 
For example, an effective momentum for a HO state can be defined by the asymptotic 
relation for large n between the radial part R„i (r) /r of the harmonic oscillator func- 
tions and the spherical Bessel functions ji{kr) of radial part of the 3D plane wave [5T| . 
Kallio showed that this relation is very accurate at small r for all n values [35]. The 
alternate definition, suggested by Vary [33] , identifies a uv regulator with the "Kallio 
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momentum" defined by this relation so that h-aitemate = y2A. This is a scale change 
only as is the definition by fiat in [34] which arrives at the same y/2 factor for their 
A. The more important distinction is the alternate definitions of the ir cutoff which 
have different functional forms. It is clear that increasing A by increasing hjj in a 
fixed- model space is not sufficient; doing so increases both of the putative infrared 
cutoffs as well because A = \\/N + 3/2 = Xsc{N + 3/2) and one continues to effec- 
tively calculate in an effective confining volume which is getting smaller rather than 
larger. This confining volume is certainly removed by letting N — >■ oo, at fixed Huj, 
because HO functions form a basis of the complete space. In addition, taking N ^ oo 
simultaneously removes the uv cutoff defined by A and the ir cutoff defined either by 
A or Asc. But increasing N without limit is computationally prohibitive. Thus there 
is a practical issue to address: whether one must take the ir cutoff to zero by taking 
huj — at fixed A (Air = A definition) or whether it is sufficient to allow Huj be some 
larger value, perhaps near that used in traditional shell-model calculations, and let 
an increasing N take Xir to small values, as it does with the definition Xir = Xsc- 

4 A study of uv and ir cutoffs in the triton 

We display in a series of figures the running of the ground-state eigenvalue of a single 
nucleus, ^H, on the truncated HO basis by holding one cutoff of (A, A^^) fixed and 
letting the other vary. These '^H calculations were made for iV < 36 and values of hut 
as appropriate for the chosen cutoff value. For N > 16, we used the code manyeff 
[25] which uses Jacobi coordinates and the No-Core Shell Model Slater Determinant 
Code [27] which use single-particle coordinates for smaller N. We checked that the 
codes gave the same eigenvalues for overlapping values of N, indicating that the 
Lawson method satisfactorily restores translational invariance to ground-state energy 
calculations in single-particle coordinates. 
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Figure 2: (Color online) Dependence of the ground-state energy of '^H (compared to 
a converged value; see text) upon the uv momentum cutoff A = ^ymN{N + 3/2)huj 
for different fixed A = \/mNhuj. The curves are fit to the calculated points. 

In Figure 2 and the following figures, \AE/E\ is defined as \{E{A,Xir) — E)/E\ 
where E refiects a consensus ground-state energy from benchmark calculations with 
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this NN potential, this nucleus, and different few-body methods. The accepted value 
for the ground state of with this potential is —7.855 MeV from a 34 channel 
Faddeev calculation [5T], —7.854 MeV from a hyperspherical harmonics expansion 
[35], and -7.85(1) from a NCSM calculation [Tj. 

For the choice of Figure 2, Xir = X — y/iriNfioj, \ AE/E\ decreases exponentially at 
fixed A, as A increases for the values of A achieved in this study. Fixed hut implies TV 
alone increases to drive A — >■ oo. Age simultaneously. The linear fit on a semi- log 
plot is extracted from the data. For fixed A, a smaller A implies a smaller \AE/E\ 
since more of the infrared region is included in the calculation. 
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Figure 3: (Color online) Dependence of the ground-state energy of "^H (compared to 
a converged value; see text) upon the ir momentum cutoff A = ^/rrifqhLj for fixed 
A = ^JmN{N + il2)nuj. 

In Figure 3 we hold fixed the uv cutoff of (A, Air) to display the running of \/S.E /E\ 
upon the suggested ir cutoff A. For fixed A, a larger A implies a smaller \/S.E/E\ since 
more of the uv region is included in the calculation. But we immediately see a 
qualitative change in the curves between the transition A = 700 MeV and A = 800 
MeV; for smaller A, \IS.E/E\ does not go to zero as the ir cutoff is lowered and more of 
the infrared region is included in the calculation. This behavior suggests that \ A.E /E\ 
does not go to zero unless A > A^^, where A^^ is some uv regulator scale of the 
NN interaction itself. From this figure one estimates A^^ ^ 800 MeV/c for the 
Idaho N'^LO interaction. 

Yet the description of this interaction in the literature says that the version we use 
has the high-momentum cutoff of the regulator set at A^vslo = 500 MeV/c \2\\ . This 
does not mean that the interaction has a sharp cutoff at exactly 500 MeV/c, since 
the terms in the Idaho N'^LO interaction are actually regulated by an exponentially 
suppressed term of the form 



exp 



In this expression, p and p' denote the magnitude of the initial and final nucleon 
momenta of this non-local potential in the center-of-mass frame and n>2. Because 
the cutoff is not sharp, it should not be surprising that one has not exhausted the 
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uv physics of this interaction for values of single-particle A somewhat greater than 
500 MeV/c. Note that this form of the regulator allows momentum transfers {p — jf) 
to achieve values in the range up to 2Ajv3lo- Can one make an estimate of the 
uv regulator scale of the Idaho N^LO interaction which is more appropriate to the 
discrete HO basis of this study? An emulation of this interaction in a harmonic 
oscillator basis uses huj = 30 MeV and Nmax = N = 20 [SS]. Nucleon-nucleon 
interactions are defined in the relative coordinates of the two-body system so one 
should calculate A^^ — ^Jm[N + 3/2)Huj with the reduced mass m rather than the 
nucleon mass appropriate for the single-particle states of the model space. Taking 
this factor into account, the successful emulation of the Idaho N'^LO interaction in a 
HO basis suggests that A^^ ~ 780 MeV/c, consistent with the figure. 



For A < A^^ there will be missing contributions of size |(A — A^^)/A''^'^| so 
"plateaus" develop as A — >■ 0, revealing this missing contribution to \AE/E\. We 
cannot rule out the possibility of a plateau appearing at the level of 0.0001 or less for 
A > 800 MeV/c as A — > 0. This is because the smallest A available to our calculations 
is hmited by A = A/^/N + 3/2 and the largest = 36 with our computer resources. 
That is, the leftmost calculated points of Figure 3 move to higher values of A as fixed A 
increases above 800 MeV/c. At fractional differences of 0.001 or less, the development 
of possible plateaus could be masked by round-off errors in the subtraction of two 
nearby numbers, each of which may have its own error. Nevertheless, the "plateaus" 
that we do see are not flat as A — >■ and, indeed, rise significantly with decreasing 
A < A^^ . This suggests that corrections are needed to A and A which are presently 
defined only to leading order in A/ A. The authors of |31] take our suggested simile 
of a truncated basis to a confining region quite seriously and use it to obtain a first 
order correction to both A and A^^. We hope to learn if higher-order corrections can 
be directly determined by our data in a future study. 
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Figure 4: (Color online) Dependence of the ground-state energy of ^H (compared to 
a converged value; see text) upon the uv momentum cutoff A = y/m]\i[N + 3/2)huj 
for different values of the ir momentum cutoff Xsc = \/ {m]sffiu})/{N + 3/2). Curves 
are not fits but simple point-to-point line segments to guide the eye. 

Now we turn to the second pair of cutoffs of (A, A^^) and display in Figure 4 the 
analogue of Figure 2 except that this time A^^ = Asc — \J rriNf^ / [N + 3/2). For 
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fixed Xsc, \AE/E\ does not go to zero with increasing A, and indeed even appears to 
rise for fixed Xsc > 35 MeV/c and A > 800 MeV/c. Sucli a plateau-like behavior was 
attributed in Figure 3 to a uv regulator scale characteristic of the NN interaction. Can 
the behavior of Figure 4 also be explained by a "missing contributions" argument; 
i.e. an argument based upon Asc < X^/^ where is a second characteristic ir 

regulator scale implicit in the NN interaction itself? One can envisage such an ir 
cutoff as related to the lowest energy configuration that the NN potential could be 
expected to describe. For example, the inverse of the np triplet scattering length 
of 5.42 fm corresponds to a low-energy cutoff of about 36 MeV/c. The previously 
mentioned emulation of the Idaho N^LO interaction in a harmonic oscillator basis |3B] 
has A^^ ~ 36 MeV/c. At low A and A.^ < A^^, \AE/E\ does fah with increasing A 
and this behavior can be fitted by a Gaussian as shown for and and other s-shell 
nuclei in [1]. But we will see in the next figure that one has not yet captured the uv 
region at these low values of A. 



0.1- 



\AE/E\ 
0.01- 



0.001 ■ 



0.0001 ■ 



[TTlMTDnD □ n □ □ □ 


□ 


• 


mccooo o o o o o 




• 




A 

• 










'^iAAi A A A a ^ " 












Idaho N^LO 








o • 

o r 






• 

O * 


A (MeV/c) 








■ 


□ 


400 


■ 


o 


500 




A 

■0 


600 


■ 


■ 


700 






800 


• 


• 


1000 


■ 

■ • 




1200 



10 20 30 40 50 60 70 80 90 

X,c (MeV/c) 



Figure 5: (Color online) Dependence of the ground-state energy of (compared to a 
converged value; see text) upon the ir momentum cutoff Xsc = \/ {mi^hu}) /{N + 3/2) 
for fixed A = ^mAr(iV 3/2)fiw. 

Figure 5 is the analogue to Figure 3: only the variable on the x-axis changes from 
A to Asc = A^/A. For A < A^^ ^ 780 MeV/c the missing contributions and resulting 
"plateaus" are as evident as in Figure 3. (Please see discussion of Figure 3 for an 
account of possible "plateaus" for larger values of A.) The tendency of these plateaus 
to rise as Asc — ^ again suggests a refinement is needed to this first-order definition 
of the cutoffs. Around A ^ 600 MeV/c and above the plot of \AE/E\ versus Asc in 
Figure 5 begins to suggest a universal pattern, especially at large Asc. For A ~ 800 
MeV/c and above the pattern defines a universal curve for all values of Asc- This is the 
region where A > A^^indicating that nearly all of the ultraviolet physics set by the 
potential has been captured. Such a universal curve suggests that Asc could be used 
for extrapolation to the ir limit, provided that A is kept large enough to capture the 
uv region of the calculation. Figure 5 is also the motivation for our appellation Asc, 
which we read as "lambda scaling", since this figure exhibits the attractive scaling 
properties of this regulator. 

We now utilize the scaling behavior displayed on Figure 5 to suggest an extrapo- 
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Figure 6: (Color online) The ground state energy of calculated at five fixed values 
of A = y/ TTiNiN + 3/2)hLu and variable Xgc = y/ {niNfito) / (N + 3/2). The curves 
are fits to the points and the functions fitted are used to extrapolate to the ir limit 
Xsc = 0. 



lation procedure which we demonstrate in Figure 6. The extrapolation is performed 
by a fit of an exponential plus a constant to each set of results at fixed A. That 
is, we fit the ground state energy with three adjustable parameters using the rela- 
tion Egs{Xsc) — aexp(— 6/Asc) + Egsi^sc = 0). The mean and standard deviation 
of the five values of Egs{Xsc = 0) were -7.8511 MeV and 0.0011 MeV, respectively, 
as suggested by Figure 7 in which the overlap of the five separate curves cannot be 
discerned. It should be noted that our five extrapolations in Figure 7 employ an ex- 
ponential function whose argument 1 /Xsc = \/ [N + 'i/2) / {rriNtujj) is proportional to 
\/N/ (hui). This extrapolation procedure of taking X^c — \/rni^fiio/ (N + 3/2) toward 
the smallest value allowed by computational limitations treats both N and huj on an 
equal basis. The exponential extrapolation in N /{huj) is therefore distinct from the 
popular extrapolation which employes an exponential in N^ax {— N for this s-shell 
case) [II HI m [20j and provides a refinement to the procedures of the 1970s for dealing 
with "smooth" potentials. 

This extrapolation procedure treats both N and hoj on an equal basis. For ex- 
ample, the extrapolation at fixed A ~ 1200 MeV/c employs values of huj from 41 to 
65 MeV and iV = 22 - 36. The one at fixed A = 800 MeV/c employs values of hv 
from 18 to 44 MeV and TV = 14 — 36. The curves of Figure 6 encompass values of 
Xsc between 20 and 52 MeV/c. We attempted to quantify the spread in extrapolated 
values by fitting only segments of the curves of this figure. Recall that the smallest 
value of Xsc requires the largest N. Fits to the segment from Xsc = 20 MeV/c to 
Xsc = 40 MeV/c (always for the five displayed values of fixed A) resulted in a mean of 
—7.8523 MeV and standard deviation of 0.0008 MeV. Cutting out the left hand parts 
of the curves and fitting only from Xsc = 30 MeV/c to Xsc = 55 MeV/c gave a mean 
of -7.8498 MeV and standard deviation of 0.0022 MeV. For both these trials a rather 
large N was needed, ranging from 14 to 36 but the extrapolation is quite stable. In 
contrast, values of Xsc higher than those shown in Figure 7, namely from Xsc — 50 
MeV/c to Xsc — 85 MeV/c, require fewer computational resources (A^ = 8 — 22). The 
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extrapolations have a mean and standard deviation of —7.792 MeV and 0.042 MeV, 
still not so far away from the accepted value of —7.85 MeV. 
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Figure 7: (Color online) The ground state energy of "^H calculated at six fixed values 
of A = ^JruNiN + 'i/2)tujj. The curves are fits to the points and the functions fitted 
are used to extrapolate to the ir limit A = y/mpffiuj = with fixed A as in Figure 6. 

Figure 3 suggests that an extrapolation to the infrared limit could equally well be 
made by taking A — )■ for a fixed large A. Instead we choose to extrapolate in huj with 
an eye to future exploitation of archival calculations made in the variables {Nmax,fi^)- 
In Figure 7 we fit the ground state energy of with three adjustable parameters 
using the relation Egs{huj) = aexp(— c/Sw) + Egs{fvjj = 0) six times, once for each 
fixed value of A. It is readily seen that one can indeed make an ir extrapolation by 
sending hw ^ Q with fixed A as first advocated in Ref. [10] and that the five ir 
extrapolations with A > A^''^ ~ 780 MeV/c are consistent. The spread in the six 
extrapolated values is about 0.049 MeV or about 1% about the mean of —7.832 MeV. 
The standard deviation is 0.020 MeV. 

Now let us accept the role of the ordered pair (A, A.^^) of cutoffs in these varia- 
tional calculations and examine the ordered pair {Af,huj). That is, we take the basis 
truncation parameter Af and the HO energy parameter huj to be variational param- 
eters. We now observe convergence as the truncation of the model space is lessened 
by increasing N = Nmax , where N is the specific truncation parameter J\f and N^ax 
is the total number of energy quanta kept in the basis. Figure 8 shows a plot of the 
variational energy of the ground state of plotted in this traditional way, pioneered 
in Figure 1 of [5] and continued through |37| to the present day [HI [31] . Optimum 
values for the parameters that enter linearly can be obtained by solving a matrix 
eigenvalue problem. But the optimum value of the nonlinear parameter must in prin- 
ciple be obtained by, for example, numerical minimization which could be difficult 
as the algorithm could easily miss the global minimum and get trapped in a local 
minima. The plots such as Figure 8 and others in the nuclear physics literature show 
that 1) for small bases a change in the non- linear parameter hco can have a dramatic 
change in the variational estimate of the ground state energy and 2) the dependence 
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Figure 8: (Color online) Dependence of the ground-state energy of "^H upon hoj = 
X^/ruN = XlJ[mNiN + 3/2)] for fixed N = A^/X^ - 3/2 = A/X,,, - 3/2. Curves are 
not fits but spline interpolations to guide the eye. 



on the nonlinear parameter decreases as the basis size increases. These observations 
seem to vitiate the need for an extensive numerical minimization by varying huj [38] . 
For example, in Figure 8 the minimum of each fixed N curve is easily read off the 
plot. 

From Figure 8, we see that the variational energy decreases and thus moves away 
from the converged value -7.85 MeV as huj at fixed N (for all N considered!). 
This is readily understood in terms of Figure 1. At fixed N one captures more infrared 
physics by lowering the infrared cutoff (A^^ oc ^/Kj) but misses the ultraviolet physics 
because lowering hu! also lowers the ultraviolet cutoff (A oc y/Huj). The loss of uv 
physics due to the lower htj overwhelms the gain of ir physics and the estimate of 
the ground state becomes very bad. A similar situation holds as hu increases: the uv 
cutoff increases toward oo so that more uv physics is captured but the ir cutoff also 
rises and more and more of the infrared physics is lost to the calculation. 

The approximate minimum of the N ~ 8 curve is at hio ^ 43 MeV which corre- 
sponds to A ~ 620 MeV/c and A^^ ^ 65 MeV/c. From Figure 5 we realize that for this 
small value of A < A^^ ~ 780 MeV/c and large value of A^^ > X^/ ~ 36 MeV/c, we 
would expect about a 30% shortfall in the ground state energy and this is what we 
see in Figure 8. At the minimum of the N — 8 curve the variational parameters are 
nowhere near their limits in the (A, A^^) regulator picture and the variational energy 
is not very good. Because N (x A^/A^ or oc A/Xsc, increasing the truncation pa- 
rameter N simultaneously increases the uv cutoff and decreases the ir cutoff so that 
the curves move lower and lower. We observe that, as fixed N increases, the minima 
of each curve moves to a lower value of ttu, as was previously observed in similar 
calculations for *He with this potential [7] and for a variety of nuclei {A = 2 — 16) 
[8] with another realistic NN potential JISP16 [39]. Maris also observes a monotonic 
movement to the left with a basis truncation on the single-particle basis so that the 
truncation parameter J\f becomes N shell rather than Nmax [10] ■ Apparently another 
behavior, first a shift to the right and then to the left as fixed N is increased, is noted 
in [31] and interpreted as first an approach to uv convergence and then, as the uv 
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physics is obtained a further convergence in the ir regulator. We, and other NCSM 
calculations (including one with a N shell truncation), do not see this behavior. 

In Figure 8, the monotonic movement to a lower is clear as N increases from 8 
to 20, all values corresponding to A < A''^^ ~ 780 MeV/c, the region in which the uv 
physics has not yet been captured. As is increased to A^ = 24 (not yet possible for 
p-shell nuclei with present day computers and codes) the minimum moves down to 
fwj ~ 24 MeV which corresponds to A ~ 790 MeV/c , A ~ 150 MeV/c and A^^ ~ 31 
MeV/c. At these values the uv cutoff seems high enough (see Figures 3 and 5) and the 
ir cutoff low enough (see Figures 2 and 4) that one could argue that convergence was 
nearly reached. As A^ increases from 24 to 36 the fixed A^ curves pile up on each other, 
but an expanded scale (not shown) separates them to demonstrate that the minimum 
stays near 24 MeV (A ^ 920 MeV/c and A^c ^ 24 MeV/c) and the curves become 
somewhat independent of huj within a limited range. Even so, any calculation in a 
finite basis should be examined from the point of view of the more physical regulators 
(A, Air)- This calculation should, in principle, always be extrapolated to the uv and 
ir limits. Independence of hu for fixed N is due to a playoff between the uv and ir 
cutoffs and it should be understood how this playoff affects the calculation. The often 
heard mantra "look for independence upon the value of hw because that means the 
series of calculations has converged" should be retired, in my opinion. 

After submission of [Ij, Furnstahl, Hagen and Papenbrock posted an investigation 
of uv and ir cutoffs in finite oscillator spaces [34] . They assume that A^c (scaled by 
a factor of from the Xsc of this paper) is the ir cutoff. They take our suggested 
simile of a truncated basis to a confining region quite seriously and use the simile to 
derive an explicit extrapolation formula in their ir cutoff. The derived formula is the 
same (exponential in N/{huj)) as the one of pj reviewed here and is used in the 
same way: establish that the uv cutoff is large enough and then extrapolate in the 
ir variable. In addition, they suggest a first (higher) order correction to both the uv 
and ir regulators. The caveat to what they call a "theoretically derived ir formula" 
is the remark made recently by Lieb et at. "If one fixes the particle number N in a 
very large box and calculates the shift in energy caused by [a given local one-body 
potential] V, the answer depends on the box shape and boundary conditions" [41 . 
But this has always been true [28]. 

As [31] assumes that (scaled) Xsc is the ir regulator, they took the behavior shown 
in Figure 4 for small Xsc to suggest a second extrapolation formula for the uv cutoff. 
That is, at low A and Xsc < X^^^ , \ AE/E\ falls with increasing A and this behavior can 
be fitted by a Gaussian, as shown for and other s-shell nuclei in [1]. This Gaussian 
in A cx yjN{fiijj) then becomes an exponential in N(huj). Their final formula assumes 
relative independence of the uv and ir extrapolations so it is a sum of exponentials with 
arguments proportional to N^huj) from the uv regulator and to -\/ N/ (hu) from their ir 
regulator. These results are a useful advance on the exponential form of convergence 
in A^ (with no mention of the role of the scale parameter huj) shown less concretely by 
the forty-year-old theorems of [3] and [1] . The authors of [31] caution, as have we, that 
results such as these should be expected only for the "smooth" potentials of [3] and [1] 
(or in their momentum space characterization: "super-Gaussian falloff in momentum 
space") such as those inspired by chiral EFT or obtained by renormalization group 
transformations. The extrapolation formulae appear to be successful in calculations 
of open shell medium to heavy nuclei (A = 74) with nuclear interactions inspired by 
chiral EFT [Hj. 

There has been a recent turn to consider other bases for expanding the trial wave 
function; bases which have a presumed better behavior at large distances than the 
HO basis which has a Gaussian falloff [13]. The most effective basis used in few- 
nucleon physics |44] . in nucleon- nucleus scattering |45j and in nuclear reactions [46j 
are the Coulomb-Sturmians. This is a complete and discrete set of the eigenfunctions 
of a Sturm-Liouville problem associated with the Coulomb potential t47j. Caprio 
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et al. have recently used this basis to make NCSM calculations of light nuclei [15] . 
They found it beneficial to link the length scale parameter bi of the Sturmian with 
the length scale b of the HO eigenfunction so as to provide a closer alignment of the 
low-n Coulomb- Sturmian basis functions with the harmonic-oscillator basis functions. 
They choose to formally truncate the Coulomb-Sturmian basis with an Nmax counting 
number. Thus they end up with the same ordered pair (A/", hui) as with the HO basis. 
However, the hu value quoted for the Coulomb-Sturmian basis is simply the huj of the 
reference oscillator length, from which the actual /-dependent length parameters bi 
are chosen to align the low-n Coulomb-Sturmian basis functions with the harmonic- 
oscillator basis functions. It therefore has no direct significance as an energy scale for 
the problem. Moreover, the Nmax truncation is difficult to interpret as an "energy 
cut" as it is for the HO basis. Caprio et al. extrapolate to an infinite basis in the 
following way: the non-linear parameter Hui is varied to obtain the minimal energy for 
the highest Af available, hut is then fixed at that value and the convergence with Af 
is assumed to be exponential (extrapolation B of [5]). This is basically the procedure 
of Delves [5], rooted in theorems of functional analysis , and is not directly related 
to the EFT inspired cutoffs discussed here. Given that neither JV nor hoj are given 
an energy interpretation in this paper, it is problematic that one can simply take 
over the arguments of [Ij or [34j to define new dimensionful uv or ir cutoffs for use in 
extrapolation. Yet the savings in computation and increase in physical understanding 
should motivate such an effort in the future. 

In summary, we have introduced a practical extrapolation procedure with A — >■ oo 
and Xir —?' which can be used when the size of the HO basis needed exceeds the 
capacity of the computer resources as it does for ^He and ^He and certainly will for 
any more massive nuclei. Unlike other extrapolation procedures the ones advocated 
in this paper treat the variational parameters Af and huj on an equal footing to extract 
the information available from sequences of calculations with model spaces described 
by {Af, huj). We have established that A does not need to be extrapolated to oo but 
if A > A^^ set by the potential one can make the second extrapolation to zero with 
either ir cutoff Xsc (see Figure 6) or A (see Figure 7). The choice of the scaling cutoff 
Xsc is especially attractive as A need not be held constant but any A large enough 
can be used in the ir extrapolation. The traditional plots in the variables (A/", hui) can 
be understood by considering the uv and ir cutoffs as primary. 
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